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A review of the scaled equation of state proposed for the critical region of fluids and magnets is 
given using the language appropriate for fluids. The experimental evidence for the validity of the basic 
hypothesis underlying this equation of state is discussed in detail. Experimental data in the critical 
regions of CO2, Xe, and He 4 are then analyzed using a closed-form expression for the chemical potential 
as a function of density and temperature, based on scaling ideas. Agreement between the proposed equa- 
tion and the experimental data is found for the three substances. The results of the scaling of A/a, Ap, t 
data are shown not to he in contradiction with the analysis, also based on scaling ideas, of independent 
experimental measurements of both specific heat and vapor pressure. 
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1 . Introduction 

The thermodynamic anomalies in the critical 
region are usually described by means of power 
laws, assumed to be asymptotically valid approaching 
the critical point. Those of direct interest to us are, in 
Fisher's notation [l] 1 : 

Coexistence curve: Ap — Bi—t)^ (1.1) 

Critical isotherm: 



\p(p, T c )-P(pc TV) I ~| M (p, T,.) -flipc, T e ) I 

= A-Ap|Ap| s -' 



(1.2) 



Compressibility: 

on the critical isochore, K T - \t~ y , T > T r 

along the coexistence curve in the one-phase region: 



p 2 Kr^T'(-t)-y T<T r 

Specific heat at constant volume: 
on the critical isochore, 



(1.3) 



C r = (A + la)t~« T>T C 

C,« (4f /«')(- *)•*' T<T < <1.4) 

along the coexistence curve in the one-phase region: 
pC-Wf/a') (-*)-"' T<T C 
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Here Ap= (p — p r .)/p t -, t= (T—T r )/T r and all proper- 
ties are measured in reduced units through the appro- 
priate combinations of critical parameters. Thus the 
temperature is measured in units T c , the density in 
units p c , the chemical potential /x(p, T) in units pdp r . 
the specific heat per unit volume pC, in units p c /T r . 
etc. The critical density, pressure and temperature 
are therefore unity bv definition. 

The scaled equation-of-state, recently proposed 
[2, 3, 4], is a partial formulation of the thermodynamic 
behavior in the critical region, incorporating these 
anomalies. In this paper we discuss the critical 
region of fluids in the light of this theory. A review of 
the basic hypothesis and the formulation of the scaling 
laws for fluids is given in section 2. The experimental 
evidence for the validity of the basic hypothesis is 
discussed in detail in section 3. Our previous tests [5] 
of the extent to which experimental data on fluids 
possess the scaling property are recalled, and the 
points where the early analysis needs improvement 
are discussed (sec. 4). An explicit expression for the 
equation-of-state in terms of scaled variables is pro- 
posed in section 5. The details of the procedure 
followed in analyzing the experimental data are given 
in section 6. The results of the analysis of the data, 
complete with internal and external consistency 
checks, are reported in sections 7 to 9. Finally, con- 
clusions are presented in section 10. 

2. Scaled Equation-of-State 

The equilibrium thermodynamic properties of a fluid 
are determined by the knowledge of either the Helm- 
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holtz free energy per unit volume a{p,T) as a function 
of density p and temperature T, or of the corresponding 
free energy per mole A(v, T) as a function of molar 
volume v and temperature. In the first case the variable 
conjugate to the density is the chemical potential 
/jl= (da/dp) t while in the second case the variable 
conjugate to the molar volume is the pressure 
P = — (dAldv)r. It is a question of which pair of conju- 
gate variables is most suitable for discussing thermody- 
namic properties in the critical region. We have made 
the choice on the basis of symmetry arguments dis- 
cussed more fully in the next section. Magnetic 
critical phenomena show, by definition, complete 
symmetry with respect to the line M = 0. Thus, the 
spontaneous magnetization is fully symmetric around 
M = 0. Since the analogy of magnetic and fluid critical 
phenomena plays such a basic role in scaling law 
theory, we choose the pair of variables which shows 
the best symmetry around the critical isochore. It 
is well known that for simple fluids the phase boundary 
shows more symmetry around p c in the pT plane than 
around v c in the vT plane. This and other symmetry 
considerations (sec. 3) lead to a preference for the 
pair p,, p. 

Widom [2], using chemical potential, density and 
temperature as variables, first proposed an equation- 
of-state for fluids; it is exactly equivalent to those 
proposed for the Ising model and for magnets [3, 4], 
as was shown by Griffiths [4]. The main assumptions on 
which this equation is based are the plausible ones that 
free energy and chemical potential are regular in p 
and t throughout the entire one-phase region with 
possible exclusion of the phase boundary. Moreover, 
the quantity Ap = p(p, t) — fJi(pc, t) is also assumed 
to be regular throughout the entire one-phase region, 
and, at least asymptotically, antisymmetric around p r . 
It should be emphasized that regularity of p and A/x 
implies regularity of pipe, t) at the critical tempera- 
ture, as is the case in the lattice gas. 

In the following, a review of the properties of the 
equation-of-state is given according to Griffiths' 
formulation [4]; then, the derived properties of the 
free energy and the specific heat are examined. 

2.1 . The Equation-of -State 

The asymptotic equation for the phase boundary, 

eq (1.1), suggests the introduction of a variable x 
defined by 

x = f/|Ap|»«». (2.1) 

In terms of this variable, it is hypothesized that the 
equation of state is of the form: 



kp = kp\kp\*- l h(x) 



(2.2) 



The quantity h(x) is a function of the variable x only, 
which according to the requirements of analyticity of 
a(p, t) , p(p, t) . and A/x(p, t) must have the following 
properties: 

(1) It is a real positive function of x everywhere in 
the range — *<><*< x which vanishes at the phase 
boundary x = — xo, xo = B~ l/l3 , with a finite slope: 



¥) = (P 2 K T )-^ -0-»A'(-*H-«-i>,«-i) 

PP /coex. 

(2.3) 

and from the definition given in eq (1.3) 

y' = /3(8-l) 

r' = [j8- 1 ft'(-*o)*if v ']" 1 - (2.4) 

About x = h(x) possesses the series expansion 

M#) = 2 h j x j =ho + h ] x + h 2 x 2 + . . . 
J=o 

so that on the critical isotherm \Ap\ = A |Ap| 6 and 
from the definition in eq (1.2) 



(2) For large x 9 h{x) has the series expansion 

h(x) = ^7]nX^ +1 - 2n} 

11=1 

which means that 



(-) 



[Kt^^V^-V 



and 



y =jB(8-l)=y' 

T^nr 1 
(3) The compressibility is given in general by 



[p«Xr]-»-|Ap| 



, 8-1 



8h{x) — - xh' (x) 



P8h{x) ^xh'(x). 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



For it to be nonnegative h(x) must satisfy the condition 



(2.10) 



Thus, with the definition *y = /3(8 — 1), the scaling 
law eq (2.2) is consistent with anomalies (1.1, 1.2, 
1.3). 

2.2. The Free Energy 

The free energy per unit volume a(p, t) is only 
partially determined by the equation of state: 
p=p{p, t)— (da/dp) t . First, since the scaling law 
equation is an expression for A/x, an unknown function 
pipe t) multiplied by p enters into the expression 
for a(p, t) . Secondly, an integration constant A (T) 
will occur additively. It is assumed that Aq(T) is 
an analytic function of the temperature. 

In the general case a#0 (0< a = 2-/3(8+l) < 1) 
we have the following expression for the free energy: 
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a(p, t) = A»(T) + pp(p<, t) + |Ap| 8+1 a a (*) 



Pv (T)=-A (T)-a a (-xo)xr 2 (-t) 2 - a ,t<0 



with a a (x) related to h(x) through: 

-xa' a {x) + (2-a)a a {x) = (3h(x) 



(2.11) 



(2.12) 



p(p c , t)=-A (T)-Ct 2 -«, t>0. 



(2.19) 



a a (x) is analytic everywhere for — jco < * < x : near 
X X Q it can be written in the form 

/ \ n\ ^0 h\X 

[2 — a 1 — a 

-*!*! 1 -* J [*(y)-Aiy-Ao]|y| a -*rfy|. (2.13) 

In the two-phase region da(x) is constant and equal 
to its value on the phase boundary a a {—xo). (The 
numerical value of a a {— xq) can be evaluated by using 
in eq (2.13) the series expansion for h(y) near y=0 
with as many terms as necessary to cover the range 
-*o,0). 



For large positive x, a a (x) can be expanded as 

1 



a a (x) = Cx 2 ~ a -\-- ^ n^rinX 2 -"- 2 * 11 

the constant C being given by 

C = -p( y [h(y)-ho-h l y]y-' i dy. 
Jo 

In the case a*0, the free energy is given by 

fl(Ap, t)-A 9 (T)+li{pc 9 t)p+|Ap|* +1 a (*) 

+ htPln |Ap| 

with ao(x) satisfying the differential equation 

— xa'o(x) +2oto(x) = p[h(x) — h 2 x 2 ]. 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



The necessary and sufficient condition for the heat 
capacity to be positive is that a" a {x) ^ for all x in the 
range — Xo < x < °°. This is hard to translate into any 
simple condition on h(x) or its derivatives. However, 
the condition that for all x in the range — xo< x < « 
h"(x) be nonnegative is sufficient for the positivity 
of the heat capacity and is useful in constructing 
functions. 

From the thermodynamic relation p = pip — a, 
through eq (2.2, 2.11), the equation of state in the 
variables p, Ap, t may be derived: 

P (p, T)=Ap\Ap\*-<h(x) + \Ap\^>[h(x)-a a (x)] 

-Ao(T). (2.18) 

The behavior of the vapor pressure p v (t) and the 
pressure above T ( on the critical isochore will then 
be given, respectively, by 



The nonanalytic term is of the form \t\ 2 ~ a . Its co- 
efficients for positive t differs from that for negative /. 
However both coefficients are related to h(x), through 
eq (2.13) and (2.15) respectively. They are both pos- 
itive since a a ( — xo) and C are negative. Since the 
chemical potential along the critical isochore is as- 
sumed non singular at t = 0, the singularity in the 
pressure is directly related to the singularity in the 
specific heat according to the scaling laws. Indeed 
the second derivative of the pressure on the critical 
isochore according to eq (2.19) diverges with the ex- 
ponent a and the same exponent is found for the 
specific heat. From 



pC r 



r& 



{T in units T c , pC v in units PrjT c ) (2.20) 

the scaling equations for the specific heat are derived 
to be: 

r IT d 2 A () , d 2 p(pc ,t) d' 2 fi(p c , t) 

+ a(-x )x« - 2 (2-a)(l-a)(-t)-« 
two-phase region (2.21) 

+ |Ap|-^a;;u) 

one-phase region. 

The singular behavior is described by the term with 
( — t)~ a in the two-phase region and by the term with 
|Ap|~ a/ in the one-phase region. It is therefore de- 
rived completely from the properties of h(x): in partic- 
ular along the critical isochore in the one-phase 
region: 



a" a {x)=C{2-a)(\-a)x-<*. 



(2.22) 



Thus pC v jT behaves as t~ a : also on any other isochore 
t~ a is the asymptotic behavior for large x. The complete 
expressions for the specific heat on the critical isochore 
are: 

C r= <PAp . (Pu(Pc J) 

T dt 2 * dt 2 
+ a a (-x Q )xS- 2 (2-a)(l-a)(-t)- a , t<0 

_Qv = d 2 A () , d 2 pi(pc ,t) 
T dt 2 ^ dt 2 

+ C(2-a)(l-o:)r«, t>0 (2.23) 
in striking analogy with (2.19). 
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The coefficients defined in eq (1.4) are found to be 
related to the properties of h(x) by: 

A + /a = -C(2-a)(l-a)=p \ '' K\y)r~ x dy 

Au/a = — {2 — a){l — a)a a {—x )x$- 2 

A l -la = pj i) h"(y)\y\^dy (224) 

while the discontinuity in the specific heat across the 
phase boundary, ApC r , is given by 



ApC r = (3xr l h'(-xv)(-t)-<*. 



(2.25) 



For the case a = 0, analogous expressions can be de- 
rived from a free energy of the form (2.16). In this case 
the specific heat has a symmetrical logarithmic sing- 
ularity on the critical isochore with a superimposed 
finite jump and the pressure at p = p c has singular 
terms of the form at 2 \n \t\ + bt 2 with the same value 
of a but different b depending on whether T^O. 
Details are given in appendix I. 



3. Scaling Symmetries and Real Fluids 

3.1. Antisymmetry of the Chemical Potential 

The scaling laws assume certain symmetries in the 
thermodynamic variables which are present in mag- 
netic systems, in the Ising model, and therefore in 
the lattice gas, but cannot a priori be expected to hold 
in real fluids. In the magnetic case, the spontaneous 
magnetization is symmetric w.r.t. the line M = 0, 
whereas the field is antisymmetric in M. For fluids on 
thermodynamic grounds alone, we have the option to 
use a p,v or a fx,p analogy to the magnetic case. The 
lattice gas has a coexistence curve symmetric in the 
density whereas the chemical potential is antisym- 
metric with respect to the critical isochore. For real 
gases, a consideration of the p,v and the fi,p behavior 
at different temperatures strongly favors the second 
option. 

In figures 1 to 3, we show the behavior of Ap = p(p,£) 
— p(p c , t) versus \v — v c \ and of | A/x| versus |Apj for 
the gases C0 2 [6], Ar [7], and He 4 [8]. Above 
the critical temperature, the chemical potential 
shows a striking antisymmetry, evidenced by near- 
coincidence of the Ap(+Ap) and Ap,(— Ap) curves 
over a range of densities up to nearly 50 percent 
from critical, and up to temperatures far above 
critical. (See figs. 1 and 2 for C0 2 and Ar.) This anti- 
symmetry is absent in the p, v representation. Plotting 
Ap versus Ap, approximate antisymmetry is observed 
in a small density range only. Thus, above TV, the 
choice of variables is unambiguously the pair p, p 
from the point of view of symmetry. It may be noted 
that the practice of analyzing PVT data for critical 
anomalies using pressure and density or volume as 
variables is still widespread; since the range of anti- 
symmetry in these variables is so limited, bias is 



readily introduced into the exponents derived in this 
way. 

Below the critical temperature, the situation is 
more involved. 

The phase boundary, along which both Ap and Ap 
are zero, is symmetric around the "rectilinear diam- 
eter" p= (p L + p G )l2. Since p¥^ p v , we choose the 
variable p rather than p c for temperatures below TV. 
With this variable, neither Ap nor Ap is antisym- 
metric over an extended range, cf. figures 1 to 3; the 
p, p representation is still the better of the two, and 
far superior to p, v representation. Of course, even 
for temperatures only slightly below TV, the entire 
range of antisymmetry observed above TV (±50% 
in density) is inside the two-phase region. 

3.2. Analyticity of fJL{p Cf t) 

If the chemical potential were fully antisymmetric, 
as it is in the lattice gas, then assuming it is regular 
everywhere except perhaps on the phase boundary, 
/jl(Pc t) has to be regular in t everywhere on the crit- 
ical isochore. This is so because, due to the anti- 
symmetry, we can write 

/4p r , t) =i[p( + Ap, t) +p(- Ap, 0], 

the quantities in brackets being chosen in the one- 
phase region and therefore being regular at all t. Since 
the antisymmetry of Ap is not complete in real gases, 
we cannot assume this regularity of /JL(pc,t) to hold 
strictly. We will examine the extent to which it is valid 
by considering the specific heat, which by thermody- 
namics can be written as 

pC v = -pT(d 2 nldT 2 ) p +T(d 2 p/dT 2 ) p . (3.1) 

In the lattice gas and in the scaling law formulation, 
where p(pe, T) is regular, T(d 2 fJL/dT 2 ) Pc gives no 
anomalous contribution to C v on the critical isochore, 
the divergence of C v being accounted for entirely 
by the divergence of T(d 2 pjdT 2 ) Pf , as discussed in 
section 2. In real systems, however, an anomaly of 
(d 2 /JL/dT 2 ) Pc may be expected [9]. In fact, even in 
the simple case of a Van der Waals fluid, 20 percent 
of the jump in C v at the critical point is caused by a 
jump in T(d 2 fJL/dT 2 ) Pc , as was shown by Barieau [10]. 

From eq (3.1) it follows that in the two-phase region 
where p and p are functions of T only, pC v is linear in 
the density, with slope — T(d 2 im/dT 2 ) and intercept 
+ T(d 2 p/dT 2 ). Alternatively, C v is linear in the volume 
with slope T(d 2 p/dT 2 ) and intercept — T{d 2 p/dT 2 ). 
Experimental evidence on the density dependence of 
C v in the two-phase region near TV is unfortunately 
scarce. The best source is formed by the C v data for 
He 4 by Moldover [11]. Analyzing his raw data we find 
that —T{d 2 fxldT 2 ) is fairly constant and equal to 
3.5±0.2 in reduced units p c /pcT c . The C v data for Ar by 
Voronel [12] and for steam by Amirchanov [13] show 
little or no dependence o{—T(d 1 pildT 1 ) on temperature 
and values of about 17.5 and 46, respectively, for this 
quantity (fig. 4). Thus there is no experimental evidence 
for rapid change or divergence of T(d 2 /xldT 2 ) in the 
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FIGURE 1. 7%c chemical potential difference A/li as a function of the density dif- 
ference Ap(a), compared with the pressure difference Ap as a function of Ap(b), 
and with Ap as a function of the volume difference Av(c), for three isotherms of 
C0 2 (ref. [6]). 

Open symbols and dashed lines correspond to the low density side of the isotherms; closed symbols and solid 
lines to the high density side. 
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Figure 2. The chemical potential difference A/x as a function of the density dif- 
ference Ap(a), compared with the pressure difference Ap as a function of Ap(b), 
and with Ap as a function of the volume difference Av(c), for two isotherms of 
Ar (ref \7\). 

Open symbols and dashed lines correspond to the low density side of the isotherms: closed symbols and solid 
lines to the high density side. 
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Fl(;URE 3. r/te chemical potential difference A/x as a function of the density dif- 
ference Ap(a), compared with the pressure difference Ap as a function of Ap(b), 
and with Ap as a function of the volume difference Av(c), for two isotherms of 
He 4 (ref [8\). 

Open symbols and dashed lines correspond to the low density side of the isotherms; closed symbols and solid 
lines to the hi§ih density side. 



two-phase region. In all these cases, on the contrary, 
rapid increase of T(d 2 p/dT 2 ) near T c is clearly 
indicated (fig. 4). 

In the one-phase region, the quantities (d 2 p/dT 2 ) p 
and (d 2 ju7dr 2 ) p in eq (3.1) are functions of density and 
temperature, thus this equation does not tell us any- 
thing about T(d 2 fJLldT 2 ) Pc . The scaling laws, however, 
predict that also in the one-phase region there is only 
one contribution to pC r which is antisymmetric in the 
density (cf. eq (2.21)). This term is linear in p and its 
coefficient is —T(d 2 jJLldT 2 ) Pc . In the cases of He 4 , Ar 
and steam [13], this linearity of the antisymmetric part 
of pC v seems to hold rather well in the one-phase 
region, and its slope varies little with temperature and 
is close to the value observed in the two-phase 
region. Specifically, for He 4 we find 3.5 for 

-T(d 2 nldT 2 ) Pc 
above T c from the linear antisymmetric part of 
Moldover's supercritical C v data, and 3.9 ± .4 from Hill 
and Lounasmaa's C v data somewhat further away from 
critical [14]. We conclude that the scant experimental 
evidence presently available suggests that (d 2 p/dT 2 ) (Jc 
depends on temperature much more strongly than 
(d 2 fildT 2 ) Pc . Thus, the hypothesis that /ut(p c , T) is 
regular at T c is not contradicted by experimental 
evidence; if very careful measurements of the specific 
heat near the critical point would show an anomaly 
in (d 2 /x/dr 2 ) Pr , then the basic hypothesis of the scaling 
laws, that both p and A^l are analytic throughout the 
one-phase region, would be invalid and the scaling in 
terms of A/x, Ap, t would only be an approximation. 




STEAM 
T r -- 374 °C 



Figure 4A. (pC v /T) n for steam (ref. [13]) is reduced units as a 
function of the reduced density differences Ap. 

The slope of the specific heat isotherms is (d 2 fJildT 2 ) u and does not show any appreciable 
variations. The intercept of the isotherms at Ap = — 1 equals (d 2 pldT 2 ) u and increases 
by over a factor 2 from 300 to 373 °C. 
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FIGURE 4B. (pC v /T) n for Argon (re/. [12]) in reduced units 
as a function of the reduced density difference Ap in the 
two-phase region. 

The slope of the specific heat isotherms is — (d-fildT-) u and does not show appre- 
ciate variations. The intercept at Ap = — 1 is {d*pfd7*)u and increases considerably. 
The isochore at Ap= l.(M)6(p = 0.533 g/cm 3 ) is thermodynamically inconsistent with 
the others. 



4. A Preliminary Check of the Scaling Property 

The proposed equation of state, eq (2.2), implies an 
interesting scaling property which can be formulated 
in a variety of ways. For instance, if the chemical 
potential difference | A/x | is divided by | t |^\ eq (2.2) 
predicts that the result is a universal function of the 
variable x only. 



where 



|A|t|/M*=G(y) 
y=|*|-0=|Ap|/|*|* 



(4.1) 



Using this representation we have shown [5] that this 
scaling property holds within certain bounds in the 
critical region of fluids. The scaled chemical potential 
differences of C0 2 [6], Xe [15], SF 6 [16], N 2 [17], and 
CCIF3 [17], obtained partially from classical PVT data, 
partially from optical density versus height measure- 
ments, fall on two branches of a universal curve. The 
representation (4.1) transforms the critical isotherm to 
infinity. This is a drawback which can be avoided by 
considering other representations. In the present 
paper we will instead scale |A/x| by |Ap| 6 , thus 
obtaining the universal function h(x). With this repre- 
sentation the critical isotherm (x = 0) is simply a point 
inside the range of definition of x. In the previous 
analysis the scatter around the curve was appreciable 
for several reasons. First of all, one should expect 
deviations from corresponding states between these 
various substances. Secondly the choice of parameters 
(3 and 8 was not optimized and it was not investigated 
whether these exponents are indeed the same for all 
substances. Thirdly, it was not certain that for all the 
points correction terms to the asymptotic form of the 
equation of state at the critical point gave a negligible 
contribution. In the case of the van der Waals gas and 
in particular cases in the Ising model, the form of the 
correction terms and the range in which their contribu- 



tion is negligible are known; however, nothing is known 
a priori for real gases. Finally, the possibility of a 
systematic experimental error in part of the data 
analyzed in reference [5] is not excluded. In this paper 
we report the results of a detailed study of the experi- 
mental behavior of fluids in the critical region in which 
we have given attention to all these questions. First 
we have not assumed corresponding states but rather 
studied the gases substance by substance. This 
severely limits the choices of gases: of the substances 
analyzed in reference [5] only for C0 2 and Xe the data 
are sufficient for a detailed analysis, to which we 
can add the recent data on He 4 by Roach [8]. 

By proposing, furthermore, an explicit expression 
for h(x) it was possible to transform the scaling 
problem into a linear least-squares problem. A good- 
ness of fit criterion could be established and the scaling 
optimized by varying the parameters 3 and T ( . Finally, 
for some of the substances, a tentative determination 
of the extension of the range of validity of the asymp- 
totic form of the equation of state was made. 

5. Proposed Form for h(x) 

The series expansion for h(x) valid for large x, 
eq (2.6), contains terms of the form x y ~ 2(in . For real 
fluids the value of the exponents y and (3 are close to 
4/3 and 1/3 respectively. A closed form for h(x), 
having all the required properties, will not be found 
in a simple real expression like the ratio of two poly- 
nomials [4J. The requirement that h(x) be zero at 
x = — Xo, will be satisfied if, following the suggestion 
of Widom [2], we write 



= ^±^! cJ) ( JC ) 



Xu 



(5.1) 



The behavior of the experimental data, using reason- 
able values for (3 and 6, shows a linear dependence of 
h(x) on (x + Xo) in an appreciable range of x including 
x = 0, suggesting that 4>(x) may be considered as a 
correction term important for large values of x. Assume: 



<P(x)=E l 



1 + E-, 



' Xo 



2/3 ~| (S — D/2-1/2/3 



Xo 



n 



= E, 



1+E-. 



X 1 Xo 

X{) 



20 



(7-0/2)3 



(5.2) 



with Eu Ez constants. The corresponding expression 
for h(x): 

m * )= M~)l 1+ M^t) J (5 - 3) 

satisfies the following requirements: 

(a) is analytical everywhere for — Xo < x < 00, 

(b) /38h(x) >xh'(x): the isothermal compressibility 
is not negative, 

(c) h" (x) ^ 0: sufficient for the positivity of the heat 
capacity, 

(d) at x = — %<), h(x) is nonanalytical but h' {— xv)=E\ 
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is finite: the compressibility on the phase 
boundary is finite and behaves like (— t)~ y , 
(e) on the critical isochore (x^ 00 ), h(x) has the 
series expansion: 



h(x) = x' ] 



V anPC~ m + V b n X~ n 



(5.4) 



The first term of this series expansion leads to the 
correct behavior of the isothermal compressibility. 
The second series on the right-hand side of this equa- 
tion, however, gives spurious terms. Since experi- 
mentally 2/3 < 1 but 4/3 > 1, the leading extra term 
x y ~ x is then the third in the large # expansion of h(x). 
The equivalent expansion of A/x in terms of Ap at 
constant t > is 



A/X=Ap 



2/»(0Ap*»+2*,(0|A P | 



»/i8 



(5.6) 



The first two terms fo(t)kp and f\(t)Ap 3 are of the 
correct form but the third term gi(t) | Ap | ll(i • Ap will 
lead to a divergence of the fourth derivative of A/x with 
respect to Ap at Ap = 0. Thus the function (5.3) fulfills 
most of the conditions imposed upon h(x) but does not 
give the correct behavior for density derivatives of 
A/x higher than the third on the critical isochore. 
However we do not consider this a serious drawback. 
Experimental A/x data near the critical isochore are 
usually beset with so many uncertainties that deter- 
mination of these higher derivatives is impossible 
anyway. 

Expression (5.3) is a form for h(x) with a small num- 
ber of adjustable parameters, namely, Xo, /3, 8, E x and 
£2, which is suitable for a linear least-squares analysis 
of the experimental data. Defining 



g(x) : 



xoh(x) 



X + Xq 



2j8/(y-l) 

= [<!>(*)] 2#(y-D (5.7) 



the assumption (5.2) for <&(x) implies linear dependence 
of g(x) on [(x + X())lxo] 2(3 according to 



g (x)=Epi<y- 



1+Ei 



x + xtA 2 ? 



Xi) 



(5.8) 



Explicit expression for the parameters A, T, F ' and the 
jump in specific heat across the phase boundary, in 
terms of the parameters which appear in the suggested 
form for h(x), eq (5.3), can be derived Trom eq (2.4) to 
(2.8 and 2.25) and are as follows; 

k=E 1 [l+E 2 ] {y - im 

r- l =xo y Ei[E t Yy- l) i# 
(T')- i ^(iifi)xB y E u rir^fi[Et']< l 'y^ 



MpC v )=M- 2 Ei(-ty 



(5.9) 



The values of the parameters A + , Aft and Af can be 
evaluated numerically from (2.24) with the explicit 
form of h(x). (See appendix II for details.) 



6. Analysis of the Experimental Data 

6.1. Evaluation of the Difference in Chemical Potential 
A/x from pv Isotherms 

The difference in chemical potential 

A/x = /x(p, t) —fi(pc, t) 

along isotherms for C0 2 [6], Xe [15], and He 4 [8] is 
obtained by graphical integration of the pressure 
volume experimental isotherms according to 



Cn(p,t) fp(p,t) 

A/x = I dfi= I vdp. 



(6.1) 



In performing the integration the knowledge of the 
lower limits, the parameters p c and p(p c , t), is crucial. 
Experimentally, these parameters are known only to 
some degree of approximation and their uncertainty 
could introduce a systematic error in the analysis of 
the data. However, since A/x(p, t) is with good ac- 
curacy an antisymmetric function of Ap this property 
can be used as a criterion to establish the position of 
the point [p c , p(pc, t)~\ along an isotherm. That is, we 
will define p c and p(p c , i) as the coordinates of the 
point with respect to which A/x is an antisymmetric 
function of Ap in an appreciable range of Ap. Of course 
the values of p c must be the same for all isotherms. At 
the same time, when the data extend over a large range 
of density, an estimate of the extension of the range of 
antisymmetry can be given. Actually this has been 
possible only for CO2 and He 4 ; the range of anti- 
symmetry extends to about 50 percent in density for 
CO2 and to about 35 percent for He 4 for supercritical 
isotherms while it appears to be smaller in both cases 
for T<T C [for T<T C A,x=^ P ,t)-^(p mex , t)=fi(p,t) 
— xx(p c , t)~\. These are rough estimates which obvi- 
ously depend on the experimental error of the data. 
Since our aim is to inspect the validity of an asymptotic 
form of the equation-of-state A/x(Ap, t) which is 
antisymmetric in Ap, we will restrict the analysis 
to points which satisfy this symmetry requirement 
within experimental error. 



6.2. Gravity Effect 

For all three substances the experimental measure- 
ments have been carried out in cells of finite height h. 
When the critical point is approached, there is a large 
density gradient in the cell; consequently, the average 
density will be different from the density at the level 
where the pressure is measured. To correct for this 
error, an a priori knowledge of the density profile, 
or equivalently, of the chemical potential as a function 
of density, would be needed. Since this is the quantity 
to be determined in the first place, this knowledge is 
not available. Thus the option left is to omit those data 
that are greatly beset with this error. The effect is 
most pronounced very near the critical point, and, in 
terms of pressure, the correction is at most p c gh. If the 
difference between the measured pressure and the 
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TABLE 1. COj, chemical potential difference Ap = p(p, t) — p(p c , t) at various densities Ap= (p — p c )lp c and temperatures T for C0 2 

The Ap. Ap data are in reduced units, the factors used for reducing I lie experimental data of ref. [6| are p,= 236.7 Amagat and pdpr — 0.3075 atm Am. ' 



T°C 


10 2 Ap 


10 4 Ap 


T °C 


10 2 Ap 


W Ap 


r°c 


10 2 Ap 


10» Ap 


29.929 


28.34 


8.2 


31.185 


4.76 

9.48 

15.03 

20.93 

-3.73 

-6.99 

- 10.57 

-15.18 

-21.00 


1.3 

3.5 

10.8 

31.0 

-0.8 

-2.4 

-4.2 

-10.4 

-30.9 


32.054 


3.86 

13.98 

22.32 

-5.11 

-11.94 

-19.05 

-24.65 

-29.21 


6.2 

34.1 

87.5 

-8.8 

-24.0 

-59.2 

-112.2 

-185.0 


30.409 


23.82 


6.0 


31.320 


4.72 


1.7 


34.720 


9.76 


81.6 




-24.22 


-7.7 




9.61 

15.13 

20.96 

-3.72 

-6.93 

- 10.58 

-15.23 

-20.97 


5.3 
15.0 

37.7 

-1.7 

-3.6 

-6.7 

-15.1 

-37.7 




21.41 
- 10.96 
- 19.95 

-27.37 


235.1 

94.0 

-207.2 

-364.2 


31.013 


13.02 


3.1 


31.523 


9.50 


8.3 


40.087 


.95 


24.6 




15.20 


6.1 




14.67 


19.1 




11.34 


276.1 




21.21 


24.2 




20.28 


42.9 




-6.98 


-157.4 




27.14 


70.6 




-8.60 


-7.6 










- 12.88 


-2.8 




- 15.02 


-20.1 










- 15.50 


-6.1 




-21.17 


-49.4 










- 18.59 


-13.9 




-26.45 


-101.5 










-21.79 


-27.3 















pressure on the eritical isochore is less than three Table 2. Xe, chemical potential difference Ap = p(p, t) — p(p c , 
times this amount, we have rejected the point. at various densities Ap = (p - p. )/p. and temperatures 1 for Xe 

The Ap, Ap data are in reduced units, the factors used lor reducing the experimental 
data of ref. [5] are p r = 1.110 g/cm' 1 and p c lp c = 51.84 (em'/g I atm. 

7. The Data 

A critical examination of the internal consistency of 
each set of data has been performed before proceeding 
to the analysis. This led to the following conclusions. 

(a) CO2 data [6]. The internal consistency is good. 
Since the cell height was unspecified and varied from 
point to point, the gravity cutoff was arbitrarily taken 
at Ap = 0.006 bar. The data for Ajjl, Ap, T are listed in 
table 1. 

(b) Xe data [15]. The cell height being 1 cm, the 
gravity cutoff is at 0.003 bar. The data have been taken 
along isotherms and along isochores. As shown in 
figure 5, there is a discrepancy between the two sets 
of measurements. The reason for the discrepancy is 
not clear; however in the isochoric measurements 
check points taken at the beginning and at the end 
of each run are in agreement, while in the isothermal 
measurements no check points are available. We will 
therefore assume that the isochoric measurements are 
more reliable. We will then analyze the A/x, Ap data 
evaluated from isotherms constructed from the iso- 
choric data. This leaves us with a relatively small 
amount of data and therefore the result of the analysis 
will be of a somewhat qualitative character. The data 
are listed in table 2. 



r°c 


10 2 Ap 


10'Ap 


T°C 


10-Ap 


10'Ap 


16.790 


3.01 


0.5 


17.590 


6.14 


9.5 




6.14 


1.5 




14.57 


31.5 




14.57 


9.5 




-0.95 


-1.4 




-4.99 


-1.0 




-2.51 


-3.5 




-9.16 


-3.7 




-3.87 


-5.7 




-13.19 


-7.4 




-4.99 


-7.1 




-13.47 


-8.3 




-9.16 
-13.19 


-15.8 
-26.9 


16.990 


3.01 


1.5 


17.990 


6.14 


14.1 




6.14 


3.6 




14.57 


43.2 




14.57 


15.5 




-0.95 


-2.2 




-4.99 


-1.9 




-2.51 


-5.3 




-9.16 


-6.1 




-2.79 


-5.4 




-13.19 


-11.8 




-3.87 


-8.3 




-13.47 


-13.1 




-4.99 

-9.16 

-13.19 


-10.8 
-23.1 
-37.8 


17.190 


6.14 


5.8 


18.390 


3.01 


9.0 




14.57 


20.6 




6.14 


19.3 




-4.99 


-2.7 




-0.95 


-2.7 




-9.16 


-8.6 




-3.87 


-11.5 




-13.19 


-16.1 




-4.99 


-14.5 




-13.47 


-17.7 




-9.16 
-13.19 
-13.47 


-29.6 
-48.6 
-50.5 



571 



/o.gr/cm- 




FIGURE 5. Comparison between isothermal and isochoric 
measurements on Xe [75]. 

Full symbols and the solid lines show the isothermal measurements. Open symbols 
show the position of points obtained in isochoric measurements; for the data at T— 16.590 °C 
these are direct experimental data, for the other isotherms they are obtained by linear 
interpolation of the p, t isochores in a range of 0.01 °C. 



(c) He 4 [8J. In this experiment the dielectric constant 
of a very thin layer of He 4 (0.025 cm) between two con- 
denser plates was measured as a function of tempera- 
ture and pressure. There are three sources of sys- 
tematic error to be discussed. The first is the gravity 
effect. Although density gradients between the 
condenser plates are negligible, the pressure was 
measured at a level above the cryostat. Even if the 
head correction to the top of the cell could be ac- 
curately estimated, an unknown head approximately 
2.5 cm above the level of the condenser from the top of 
the cell, is present. This gives a gravity cutoff of 
0.0003 bar. The second source of error lies in the con- 
version of dielectric constant to density using a 
Clausius-Mosotti relation with constant molar polar- 
izability a = 0.1230 cm 3 /mol. However, if a depends on 
the density as observed for He 3 [18], 

a(cm : ymol) = 0.1234-l.ll X 10- 4 p/p r , 

then the first order correction to a density difference 
Ap calculated assuming constant a would be 
1.11 XlO 4 (l/a-fAp)Ap. This correction is different 
at the low- and at the high-density side, and would, 
at Ap ~ 30 percent, amount to 2.6 X 10 4 in Ap on the 
vapor side and 2.8 X 10 -4 on the liquid side, respec- 
tively. Since the uncertainty in Ap is about 3 X 10~ 4 , 
this is a borderline correction; however, if applied, it 
would slightly improve the antisymmetry of the 
Ap,, Ap data. 

The third source of error is more elusive and is 
related to the use of the T 58 He 4 vapor pressure scale 
[19] near the critical point of He 4 . This scale relates 
the vapor pressure to the temperature by an equation 
which is analytic at the critical point. This is contrary 
to the current ideas about critical point anomalies. 



According to the scaling laws, the divergence of the 
specific heat imples a similar divergence in d 2 p v jdT 2 . 
Moldover [11] estimates from his specific heat measure- 
ments that vapor pressures near the critical have to be 
assigned temperatures up to 0.0007 K below those on 
T 58 . One can avoid this difficulty by calibrating the 
thermometer independent of the T 58 scale, or by cali- 
brating it on this scale in a temperature range not too 
close to T c (say below 5 K). 

In the present experiment it appears that a germani- 
um thermometer has been calibrated on the T 58 scale 
in the range from 5 K to the critical point, and, more- 
over, that this calibration was extrapolated to tempera- 
tures above critical. From Moldover's estimate of the 
size of the correction below T ( , and assuming that all 
temperatures are measured relative to critical, it 
seems corrections to T c — T will be below 0.5 mK. 
For T > TV, where extrapolation is involved, the 
correction may well increase; we will assume that for 
isotherms further from critical the temperature is not 
known to better than 1 mK. 

The Apt, Ap, T data for He 4 are listed in table 3. Other 
data on He 4 are available in the work of Edwards [20]. 
These data are restricted to the region T < T c and yield 
information on the phase boundary and the behavior 
of the compressibility for t < 0. We have chosen to 
check the scaling properties on Roach's data because 
in this case also the region above the critical tempera- 
ture has been explored. We will compare the results 
with the results of Edwards only qualitatively because 
a direct comparison of the coexistence curves in the 
pT plane by the two authors shows some as yet 
unresolved discrepancies. 

8. The Analysis of the A/x, Ap, / Data 

Using conventional linear least-squares techniques, 
we fitted g(x), defined by eq (5.8), as a function of 
[(x-\- xo)Ixq] 2/3 and obtained estimates for £\, E 2 and 
their errors. The computations were performed 
using the general purpose program OMNITAB [21]. 
The experimental values of Ap,, Ap, and T (tables l - 3) 
were used to evaluate x and g(x) by using suitable 
values for p c , /3, 5, xq, and T c . We took /3 and xo from 
a preliminary analysis of coexistence curves. A more 
refined analysis, to be published [22], yields values for 
/3 and xo which are only slightly different from the ones 
we used (see table 4). The fit of the Ap,, Ap, t data is 
repeated for different values of 8 and T c until a min- 
imum in the standard deviation of the fit is reached. 
The value for T c corresponding to the minimum agrees 
with the critical temperature obtained from fitting the 
coexistence curve [221 (table 4). 

Weights were assigned to Apt, Ap, and t in the fol- 
lowing way. The precision of the experimental pres- 
sure, temperature and density was estimated from the 
description of the respective experimental methods. 
From the pressure error cr ; ,, the error in Ap, was esti- 
mated as cr l j L = 2cT l ,/pc. The estimated errors in the 
reduced quantities Ap,, Ap, t are listed in table 5. 
The fit was performed assuming that the variable 
[{x-\-xo)lxo] 2 P is free of error and that the variable g(x) 
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Table 3. He 4 , chemical potential difference Ap = p(p, t) — /*(#., t) at various densities Ap = (p — p c )lpc and temperatures TforHv 4 
The Ap, A/i data are in reduced units, the factors used for reducing the experimental data of wl\ [8| are p, =0.069 g/cm 3 and p r /p r = 2.475 10 4 torr cm'/g. 













1 1 
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TK 


10 2 Ap 


lO^A/it 


TK 


10- Ap 


KPA/a 


T K 


l()-A/> 


10 3 A/x 


TK 


10 2 Ap 


l() ! A/x 


5.1097 


33.58 


2.68 


5.1793 


-20.28 


- 1 .33 


5.1929 


13.67 


0.84 


5.221 


8.04 


1.29 




35.38 


6.71 




-21.99 


-2.29 




17.78 


2.10 




10.99 


2.30 










-25.30 


-5.08 




21.28 


4.15 




16.65 


4.66 










-28.68 


-9.74 




25.10 
28.35 


7.89 
12.91 




21.54 
25.74 


8.69 
14.38 


5.1406 


29.86 


3.77 










31.23 


18.93 




30.22 


23.47 




32.14 


7.74 










34.71 


28.95 




-8.55 


-1.40 




34.15 


12.36 


5.1795 


18.64 


0.95 




- 16.46 


-1.33 




-14.14 


-3.05 




-29.99 


-3.51 




19.64 


1.13 




-20.46 


-3.12 




-20.16 


-6.86 




-33.33 


- 10.25 




21.68 


2.40 




-24.03 


-6.11 




-25.14 


-12.88 




-36.01 


-17.57 




23.25 
31.61 


3.64 
15.63 




-28.06 
-31.77 
-35.83 


-11.93 
-20.26 
-33.83 




-30.32 


-23.71 


5.1666 


22.75 


1.25 






















25.73 


3.79 


5.1822 


19.99 


1.74 


5.201 


12.91 


1.10 










27.64 


6.14 




21.71 


2.90 




17.45 


2.62 










29.78 


9.53 




23.23 


4.09 




21.20 
25.45 


5.11 
9.76 


5.264 


1.55 

7.22 


0.69 
3.27 


5.1663 


-23.65 


-1.50 


5.1835 


-18.13 


-1.18 




29.19 


16.08 




12.90 


6.61 




-26.65 


-4.39 




-20.30 


-2.19 




33.13 


26.18 




20.14 


13.69 




-30.17 


- 10.02 










-12.03 


-0.98 




27.84 


28.12 




-33.71 


-18.75 










-15.41 


-1.98 




-3.68 


-1.55 




-35.78 


-25.61 


5.1862 


18.96 


1.86 




- 19.86 


-4.17 




-8.83 


-3.97 










21.19 


3.27 




-23.35 


-7.37 




-14.59 


-7.67 










23.65 


5.41 




-27.03 


-12.35 




-20.04 


-13.24 










-17.61 


-1.16 




-31.01 


-20.89 




-26.10 


-23.87 










-19.78 


-2.17 




-35.09 


-33.85 


5.315 


5.97 

4.54 

-9.90 

-11.12 


4.93 

3.34 

-8.04 

-9.66 



Table 4. Critical parameters for C0 2 , Xe, and He 4 . 

In most cases the error bars denote standard errors obtained by the least squares analysis. The ones marked by asterisks have been established by observing thai systematic deviations 
from the fitted equation are present if the parameter is varied by this amount or more. 





C0 2 


Xe 


He 4 


lvalues used in the fit 


|0.35 
10.135 
236.7 Amagat 
f 30.96 ± 0.04 °C* 
4.6±0.1* 
2.36 ±0.02 
I 0.30 ±0.02 
( 0.347 
0.139 
1236.8 Am 
[ 30.96 °C 


0.35 

0.186 

1.110 g/cm 3 
16.58 ± 0.03 °C* 

4.6±0.1* 

2.96 ±0.07 

0.37 ±0.03 

0.351 

0.189 

1.107 g/cm 3 
16.59 °C 


0.359 


p c from antisymmetry 


0.36 

0.0693 g/cm 3 

5. 1884 ±0.0008 K* 

4.45±0.10* 


T 1 

pr 1 

I-. >from fit of the Ap, A/x, t data 


1 1 
> from fit of coexistence curve, [22| 


2.78 ±0.03 
0.48 ±0.03 
0.355 
0.366 


Pc\ 


0.0690 g/cm 3 
5.1887 K 



contains all experimental uncertainties cr M , o" p , a t . 
Using propagation of errors the weight to be attributed 
to g(x) is found to be 



W-- 



I 



L 2/3 J iUm/ Up. 



f3(x + x„)] 



\t) \x + xj j 



[g(x)V 2 . 



\X + XoJ 

In order to obtain reliable values for the critical expo- 



Table 5. Estimated experimental errors in the reduced density 
difference (er p ), reduced chemical potential difference (cr M ) and in 
the reduced temperature (cr t ) 
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3.3X10 4 
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6.5 XIO" 5 
3.5X10 •> 
1.2X10 ' 
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3.4 x 10 ,; 


He 4 


5.0X10 5 
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nents and coefficients, we have to make sure that 
contributions from possible corrections to the asymp- 
totic form of the equation of state are negligible in 
the range of data considered. Correction terms which 
do not possess the antisymmetry property are readily 
excluded by restricting the analysis to the range of 
antisymmetry. Because the form of antisymmetric 
correction terms is unknown, their contributions can- 
not be evaluated. However, such terms will scale in 
a way different from the asymptotic term; therefore, 
points for which their contribution is important will 
show systematic deviations from the fit. For one sub- 
stance, namely C0 2 , for which many data are available 
on isotherms substantially higher than critical, the 
effect of points outside the range of validity on the 
values of the parameters may be observed and an 
estimate of the range derived. 

Actually, for C0 2 , omitting points within the range 
of antisymmetry for the isotherms between 29.929 and 
34.720 °C does not affect the results of the fit while 
including points with Ap > 11 percent at 40 °C and 
Ap > 1.4 percent at 50 °C (next point at 15%) introduces 
systematic deviations and an apparent lowering of the 
value of 8. 

For Xe the temperature range of the data is quite 
small and no appreciable effect was observed omitting 
points at large densities on any isotherms. For He 4 
omitting points within the range of antisymmetry for 
the isotherms between 5.109 K and 5.202 K does not 
affect the results while the inclusion in the analysis 
of the points with Ap > 30 percent at T= 5.221 K, 
Ap > 28 percent at 7=5.264 K and with Ap > 11 per- 
cent at T= 5.313 K introduces systematic deviations 
which again apparently lower the value of 8. 

In table 4, the values of the parameters obtained by 
fitting the A/x, Ap, t data are listed. The value of p c is 
that with respect to which Ap, is antisymmetric. The 
parameters E u E 2 are least-squares estimates, as 
explained before; the errors listed are standard de- 
viations. The parameters 8 and T c were varied step- 
wise. The values given minimize the standard deviation 
of the fit. The errors in 8 and T c are estimates obtained 
by observing that systematic deviations from the fitted 
equation arise if the parameters are varied by this 
amount or more. They do not take into account the 
effect of small variations in /3 and Xq. For the three 
substances, the values found for the exponents (3 and 
8 are the same within the error bounds. To decide if 
the small observed differences are significant, more 
accurate determinations of the exponents are needed. 
Of the other parameters describing the critical point 
behavior, E\ and Ei do not show large variations either; 
the only parameter varying greatly from substance to 
substance is Xq. 

In figures 6A to 8A the optimized function h(x) 
is compared with the experimental data for the 
three substances in a log log plot. It is seen that the 
log h(x) is a smooth and nearly linear function of 
log (x + Xo). Deviations from linearity are important 
only for quite large value of x. The point at x = (critical 
isotherm) behaves like any other point at finite x. 
The agreement with the experimental points is best 



seen from the deviations plots shown in figures 6B 
to 8B where the relative difference between the 
experimental and the calculated value of h(pc) is plotted 
versus log [(x + Xo)lxo]. No evidence of systematic 
deviations is found. 

Figure 8B indicates that for He 4 the error limits 
are somewhat too narrow. Some of the errors in table 5 
must have been underestimated. An adjustment of 
the estimated error in the pressure could remedy this 
problem. In view of the arbitrariness of this procedure, 
we have not done this. Moreover, statistical estimates 
for parameters and their standard deviations are 
usually not very sensitive to the exact weighing 
procedure used. 

From the values of /3, 8, E\ and E 2 , and using the 
equalities 



y— j3(« — 1) 

2-a=/3(8 + l) 



(8.2) 



and definitions of A, T, F' . . ., eq (5.9), we have 
determined the parameters listed in table 6. If yS and 8 
are constant from substance to substance, obviously 
the same would be true for the exponents y and a. 
Of the other parameters also A and the ratio 17 F 
seem to be remarkably steady. Therefore the critical 
isotherms of the three gases should approximately 
coincide; our previous study, implicitly assuming 
this coincidence, yielded a value of 8 slightly higher 
than the one found presently. This must be due to the 
inclusion in the previous analysis of other data not 
studied here. The compressibility ratio TIT' is sub- 
stantially smaller than our previous estimate obtained 
by graphical determination of slopes of isotherm 
(see below). Our new value agrees reasonably well 
with the value obtained for Xe by Wilcox and Balzarini 
[23] using optical interferometry, namely r/T' = 4.9. 

For each substance, internal consistency checks 
have been made by comparing experimental critical 
isotherm and compressibility data with values pre- 
dicted by our equation. 

An example is figure 9 for C0 2 , where the inverse 
compressibility (9A) on the critical isochore above T c 
and along the phase boundary below T c is shown; 
moreover, the critical isotherm is plotted (9B). For 
the latter, it is seen that the experimental data agree 
quite well with the calculated curve. This is also the 
case for the compressibility above T c . However, below 
7V, the experimental points show systematically lower 
compressibilities than our equation predicts. These 
points have been obtained by graphical determination 
of the slope of isotherms near the phase boundary, 
and by assuming linearity of p versus Ap in an appre- 
ciable range. The optical measurements [17] indicate 
strong curvature of these isotherms, suggesting that 
slopes taken from PVT data tend to underestimate the 
limiting compressibility. The function (5.3) predicts 
a very small range of linearity; the discrepancy is then 
explained by the fact that the experimental points do 
not define the true limiting behavior. The refractive 
index measurements of Edwards [20] in He 4 permit 
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Figure 6 A. Lax-log plot ofh{\) versus (x+ x )/x /br C0 2 . 

The solid line represents h{x) according to eq (5.3) with parameters given in table 4. Experimental points are 
shown with estimated errors in the variable h(x). 
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Figure 6B. Relative deviation Ah(x)/h(x)= [h^,!*) - h/i7^rf(x)]/h Pxp (x) versus 
/o#[(x + x )/x ] for C0 2 . 

Vertical error bars represent all experimental uncertainties. 
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shown with estimated errors in the variable h(x). 
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shown with estimated errors in the variable h{x). 
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another estimate of the range in which p is linear in 
Ap and also in this case important deviations from 
linearity set in at reduced densities only a few percent 
away from the phase boundary. The values of y' = 1.26 
and P = 3.8 X 10 2 determined by Edwards are in good 
agreement with our values for He 4 shown in table 6. 

Table 6. The critical exponents y and a and the calculated param- 
eters which describe the behavior along the preferred curves: 
critical isochore (F and F) and critical isotherm (A); a/50 the 
parameters describing the specific heat along the critical isochore 
[A\Afi,ApC ¥ (-t)«l 





C0 2 


Xe 


He 4 




1.26±0.02 
0.04 

5.26 X10- 2 
1.19X10- 2 
4.4 
2.6 
4.24 
5.80 
39 


1.26 ±0.06 
0.04 

5.9X~10- 2 
1.43 X10" 2 
4.1 
3.3 
3.70 
4.80 
28 


1.24 ±.04 




0.05 


r 


1.30X10" 1 


r 


3.59 X10 2 


r/r 


3.6 


A 


3.2 


A+ 


1.06 


A~ 


1.42 


^11 

AoC, • ( — t) a 


7.3 







We note that the values 0.76 and 0.56 reported tor A in ref. [5] actually were the values 
of log A. We note also that the values for ApC,( -t) a reported in ref. [28] correspond to the 
actual values multiplied by a. 
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Figure 9. Comparison between the experimental results and the 
prediction of the fit for the compressibility on the critical isochore 
above T c , the compressibility on the phase boundary below T c 
(9A) and the critical isotherm (9B). 

The experimental results for the compressibility are obtained assuming that the a, p 
isotherms are linear in p in an appreciable range around p r {T > T r ) or near the phase 
boundary (T < T e ). The solid lines represent the prediction from the fit of the Ap, Ap. 
t data. 



8. Consistency Checks: Vapor Pressure and 
Specific Heat 

In section 2 it was shown that the scaled equation 
of state predicts the nonanalytic behavior of the free 
energy a(p, t) and therefore of the pressure and of 



the specific heat at constant volume. Thus external 
consistency checks of the scaling of the A/z, Ap, £data 
are obtained by comparing the anomalies implied for 
pressure and specific heat with an independent 
analysis of the directly measured values of p and C v . 

8.1. Vapor Pressure and Pressure Along the Critical 
Isochore at T > T c 

The prediction of the scaling law for these quantities 
has been given in eq (2.19). If the analytic part^ (r) 
were known, a direct analysis of the pressure would 
yield information on the terms with | t\ 2 ~ a and on the 
value of the exponent a itself. In a relatively small 
range of t we can assume that a polynomial in t gives 
an accurate representation of the function Aq(T). 
With a given set of experimental data one can then 
observe the difference in the fit obtained with a 
polynomial in t of degree n (perfectly analytic behavior) 
or with a polynomial of degree (n — 1) plus a term 
D\t\ 2 ' a . However, the indications are that a is small 
(— 0.1); thus unless very accurate determinations of 
the vapor pressure over a not too limited range of 
temperatures are available, the values of a and D 
will have large uncertainty, even if the nonanalytic 
expression gives an improved fit to the data. The 
product Da is expected to be approximately constant 
and perhaps somewhat better determined. Vapor 
pressure data and pressure on the critical isochore 
above T c should be independently fitted and compared 
to see whether the analytic part and the values of a are 
the same. For Xe a reasonable amount of information 
is available at T <T C in the data of Michels et al. 
[24], namely, eleven data points between 247.2 K and 
289.2 K. We have fitted these data with a 3d degree 
polynomial and with the nonanalytic expression 



P-Pc 



= at+bt 2 -a a (-Xo)x%- 2 (-t) 



(9.1) 



Although the two functions have the same number of 
adjustable constants, the nonanalytic expression 
has a standard deviation of a factor 4 smaller for 
any a between 0.04 and 0.08. Due to the fact that 
there are no points very close to T c , small variations 
in T c (and p c ) do not affect the results. Values of the 
constants a, b, a n {—Xo)x$~ 2 are given in table 7. 
Because of the scarcity of data we did not study the 
dependence of the parameters on the range. 



TABLE 7. Constants of the vapor pressure equation: 

Ap/p c = at + bt 2 -a a (-xo)x«- 2 (-t) 2 - Q , 

for Xe ref [24] and C0 2 ref [25] 



a 


a 


-b 


— a a {— Xo)-x$ 2 


— aa a ( — Xo)xd~ 2 


Xe 0.04 
0.06 
0.08 

C0 2 0.04 
0.06 
0.08 


6.02 ±.01 

6.02 

6.02 

6.980 ±.002 

6.984 

6.988 


76±2 

47 
32 

85.1 ±.8 
51.1±.4 

34.2 ±.3 


82±2 

53 

38 

93.8 ±.7 
60.1 ±.5 
43.4 ±.4 


3.28 ±0.08 

3.18 

3.04 

3.75±.3 
3.60 ±.3 
3.46 ±.3 
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Table 8. Parameters describing the behavior of the specific heat along the critical isochore: cowoarison between values derived {a) from fit 
of the Ap, A/u, t data, (b) from fit of vapor pressure data, {c)from direct measurement o/C v , (d) our analysis of the experimental C v data 
ofref.lll] 



In the table blank spaces indicate that a numerical value for the pa 
on the parameter cannot be derived from the corresponding analysis. 



r could be derived by the measurement but is not available at present. A dash ( — ) indicates that informant 
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(a) 0.04 
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(a) 0.04 

(b) .04 
.06 
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2.9 
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-44 
-20to-6 
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170 
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94 

23.5 to 9.5 


39 


Xe 


28 


He 4 


7.3 
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For CO2, precise vapor pressure data have been 
recently measured [25], with 29 points between — 5 °C 
and 31 °C. Here, the data were fitted to a 3d degree 
polynomial and a nonanalytic expression analogous 
to (9.1). Use of the nonanalytic expression decreases 
the standard deviation of the fit by a factor 7. T ( , oc 
and the range were varied and an optimum fit was 
obtained for a between 0.04 and 0.10 for the range 
4 to 31 °C. From the coefficient of the nonanalytic 
term, and using the scaling laws eq (2.19) and (2.24), 
values for A^/a are derived and compared with those 
obtained from scaling the A/x, Ap isotherms of Xe 
and C() 2 in table 8. In view of the large uncertainty 
in the value of a and the correlation between a and 
A^ we feel the agreement is reasonable. In the same 
table, values for A^/a from direct measurements of 
the specific heat are also presented. We will discuss 
the agreement below. 

The constant b in the vapor pressure equation yields 
an estimate of — Td 2 Ao(t)/dt 2 , which is also given 
in table 8. 

8.2. Specific Heat 

From the scaling law eq (2.21) it is seen that in the 
one phase region 



P C V /T+ 



(PApjt) L d 2 n( P c,t) 



dt 2 



dt 2 



] iAph« 



-al(x) (9.2) 



that is, the quantity on the left-hand side is a function 
of the variable x only; its properties are related to 
the properties of h(x). If functional forms of a" a {x) and 

—7-^ could be postulated this relation could be used 

to check the scaling of data measured at various 
densities and temperatures. The situation would then 
be analogous to the scaling of the Ap,, Ap data, but 
while in the latter case only one function h(x) had to 



be chosen 

d 2 A {) 
dt 2 



ire we have two unknown functions 

and a'a(x). In order to avoid making additional 

assumptions, we decided to analyze only those specific 
heat data that are close to the critical isochore. For 
these the scaling laws predict the following behavior: 



C v =C +(A+la)Tt— 

C v ^C +(Aula)T\-t\ 
with 



t > One-phase region 

(9.3) 
t < Two-phase region 



Co = -T 



d 2 fi(p c ,t) _^_d 2 A () (t) 



dt 2 



dt 2 



assumed to be constant in a small range of t. If this 
asymptotic form is assumed valid over a comparable 
range oft above and below T c , one can determine the 
values of the parameters a, Co, A + , A~ x by fitting 
expressions (9.3) to the experimental data; systematic 
deviations should appear at the same values of \t\ on 
both branches. 

Experimental data are available for He 4 in the 
measurements of Moldover [11] and for Xe in the 
measurements of Edwards, Lipa and Buckingham [26). 
We made an analysis of the raw data on helium and 
obtained the range of values given in table 8. As an 
example, a plot of (C v - Co) T~ « with C = - 0.7 J/cm 3 K 
or 15.5 reduced units is given in figure 10. The data at 
T <T C are on a straight line of slope a = 0.06 and the 
same value of a is compatible with the data at T > T c ; 
in the latter case the points with t < 10 s may be 
slightly affected by gravity. From the values of Co and 
the known value of Td 2 fi(p c , t)/dt 2 values for d 2 Aoldt 2 
are derived. (Table 8.) The C v data for Xe have been 
analyzed by the authors [26] according to eq (9.3) with 
T assumed constant and equal to T c . They conclude 
that any value of a between and .12 gives a reason- 
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FIGURE 10. Log-log plot of the experimental data on the specific 
heatofUe\ref. [11]. 

In this plot C v . Co are in J/cm :i K. Co has been assumed constant and equal to — 0.7 J/cm :! 
K. and T c = 5.1889 K. The solid lines are lines of slope a = a' = 0.06. Similar plots are 
obtained with C„ in the range —0.3 to —0.9 J/cm 3 K. The corresponding values of a will 
range from a = 0.04 to a = 0.11. 



able fit to the data; their values for the parameters are 
also listed in table 8. 

The comparison with the values of the same param- 
eters as obtained from the analysis of the Ap, Ap, t 
data is quite good. It should however be mentioned 
that in the case of helium our value of a does not seem 
to agree with the value given recently by Moldover by 
a scaling-law analysis of his data. In appendix I we 
discuss Moldover's procedure and the possible reason 
for the disagreement. 

In conclusion, there are presently no indications of 
serious inconsistencies between the scaling-law 
approach of the equation of state, and the behavior of 
vapor pressure and specific heat of CO2, Xe and He 4 . 

9. Conclusions 

In applying the scaling laws in their present formula- 
tion to the critical region of fluids, the antisymmetry 
of the chemical potential with respect to the critical 
density is a prerequisite. We observed this antisym- 
metry in the experimental data in the critical region 
of a number of fluids in an appreciable range of density 
and temperature. As we have pointed out (sec. 3), a 
corollary of the antisymmetry is the analyticity of 
/jl(Pc T) at T c . This analyticity has been found not to 
be contradicted by the scarce experimental C v data 
presently available. Scaling of the experimental Ap, 
Ap, t data for He 4 , C0 2 and Xe has been shown to be 
possible within experimental precision in a density 
range of about ±30 percent from critical and for tem- 
peratures within — 1 percent to + 3 percent from criti- 
cal. A form we propose for the scaling function h(x) 
fits the data for these three fluids and predicts reason- 
able values for the anomalous part of C v and of the 
vapor pressure. The values of the critical exponents 
determined by the use of the function vary only slightly 
from substance to substance. There is hardly if any 



dependence on the quantum parameter A* [27]. The 
coefficient E\ and E 2 in our equation do not vary ap- 
preciably either from substance to substance; thus the 
deviations from corresponding states are almost com- 
pletely described by the parameter Xo, or B, alone, 
which parameter varies considerably from substance 
to substance. 

The present analysis has made it clear that new 
experiments are needed; in fact, several of the above 
conclusions are tentative because of lack of sufficiently 
precise and extensive experimental information. In 
particular, the fundamental question of the antisym- 
metry of Ap, should be further investigated. Most useful 
would be studies of the density profile in the gravity 
field both above and below T c , since this is a direct 
probe of the p, p relation. Particular care should be 
given to avoid vertical temperature gradients. At the 
same time, the behavior of d 2 fi(p c , t)/dt' 2 should be 
elucidated by precise C v versus p measurements in 
the two-phase region approaching the critical tem- 
perature as closely as possible. C v measurements in 
the one-phase region, particularly at T c and along the 
two-phase boundary below T c , will, in addition, provide 
better information about the exponents a and a . The 
specific heat measurements should be supplemented 
with measurements of the pressure on the critical 
isochore both above and below T c to check whether 
the divergence of C v is indeed fully explained by the 
divergence of (d 2 p/dt 2 ). Finally, precise PVT data are 
needed, in particular near the two-phase boundary. 
The most useful approach would be to obtain as many 
of these properties as possible for the same substance. 

At the present level of experimental precision the 
scaling idea, including the assumption of antisymmetry 
of Ap, and the equation we propose, seem to hold 
satisfactorily. If, from more refined data, systematic 
deviations from our function were found but scaling 
is still observed, other closed-form expressions for the 
function h(x) could be considered. For instance, a 
modification of our function, which could give a better 
agreement close to the phase boundary, and at the 
same time move the nonanalyticity at x = — xo to a 
point x = — X\(xi > xo) inside the two-phase region, is 

readily obtained by replacing!- — — J in eq (5.3) by 
J . More precise data could also reveal devia- 
tions from scaling in the present density-temperature 
range. This might lead to a refinement of the defini- 
tion of the asymptotic range without necessitating 
fundamental changes in the scaling ideas. 

However, should it be found that the property of 
antisymmetry is not fulfilled, then basic changes in the 
scaling formulation are needed, for instance by choos- 
ing new variables such as a combination of pressure 
and chemical potential. Within the framework of the 
scaling laws as presently used, the most significant 
need, from a phenomenological point of view, is the 
formulation of the free energy as a function of density 
and temperature, since the knowledge of the Ap,, Ap, 
t equation of state by itself yields only a partial descrip- 
tion of the critical region. 
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10. Appendix I 

When a = 0, eq (2.16) and (2.17) must be used to 
derive expressions for the equation of state in terms 
of the variables p, p, t and for the specific heat. 

For the pressure one obtains 

p = Ap|Ap| 8 - 1 /iU) + |Ap| 8+1 [/i(x)-a (x)] 

-h 2 t 2 \n\Ap\-A {) (T) (AI.l) 

which gives on the critical isochore: 

p(p c , t)=-A () (T)-(3Gt 2 -ph 2 t 2 \nt 

p v (T) = — A () (T) — t 2 [a (— x )xo 2 — /3/i 2 In *o] 

-ph 2 t 2 In |*| (AI.2) 

where G is an undetermined constant and a{— Xo) 
may be evaluated from the expression for ao(x) valid 
near x — 0: 

ao(x) = /3{-/io-f h\X-\-Fx 2 



\h(y) - ho- h x y- h 2 y 2 ]\y\- 3 dy} (AI.3) 
where F is a constant related to G by 



- C+J r*t 



h(y) —h () — h\y h 2 



y 



\+y _' 



(AI.4) 



The nonanalytic part is then composed of two parts, 
at 2 In t term with the same coefficient above and below 
7V, and a t 2 term with different coefficient for t > 
or t < 0. 

For the specific heat one obtains for t > 



pC v d 2 A {) d 2 p(p r , t) 
P- 



dT 1 



di l 



-ag(*)+2A 2 ln|Ap|. 

(AI.5) 



On the critical isochore: 

ao'0O=2/3A 2 ln* + 3/3A 2 +2j3G 



(AI.6) 
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in the two-phase region: 










pC v 
T ~ 


(PAo d 2 

Z dF+ p 


V>(pc, 

dt 2 


— + 2/3/i 2 


In |*| 








-2/3/* L 


In Xo - 


-f 3/3/i2 + 2 




*()) 


(AI.8) 



Therefore the asymptotic expressions for the specific 
heat on the critical isochore are 



C r =Alnt + B + t>0 



(AI.9) 



C v =A]nt + B- 



t<0 



tlth 



A = 2/3h 2 =ph"(0) 

B + = 3/3/i 2 + 2/3G (ALIO) 

B_ = 3/3h 2 - 2/3/i 2 In x + 2a (- x )x ' 2 

Except for B + , which contains the undetermined 
constant G, these coefficients can be numerically 
evaluated from the properties of h(x). 

The expression for the jump in the specific heat 
across the phase boundary is 



ApC r = (#*o)A'(-* ). 



(Ai.ii) 



In contrast to the case a > or a < 0, this disconti- 
nuity is independent of the temperature. 

Since the equation for the discontinuity in pC , 
when crossing the phase boundary does not involve 
any unknown constant, the experimental determina- 
tion of this jump would provide the most direct 
information about a . 

11. Appendix II 

Numerical evaluation of the coefficient for the 
specific heat. 
A + /a— This coefficient is given by the integral 



I> 



(y)Y* 1 dy 



(AIM) 



where h"(y)y a ~ 1 diverges neary=0. 

Assuming e small but different from the integral 
between € and infinity is evaluated by standard 
numerical integration. For the interval — e, h" (y) 
is given by 



h"(y) = 2h 2 + 6h 3 y+ 12h 4 y 2 + , 



(AII.2) 



th( 



fe € a e 

h"(y)r- x dy = 2h 2 — + 6h., 

J<> OL 



OJ+1 



+ 12A.4 



a + 2 



+ , 



(AII.3) 



For e sufficiently small, € ^ 0.Lr o , evaluating the first 
two terms only is sufficient. 

Aula — The calculation of this coefficients involves 
the evaluation of a a (— *o), which from 2.13 depends 
on the integral 



I 



(h 2 y 2 + hzy 3 + t U y 4 +. . .)\y\ a -*d\y\ (AII.4) 

{h->-h-,\y\+h i \y\ i +. .) |y|«-'d|y| (AI1.5) 
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which is of the same form as the integral involved 
in the evaluation of A + la. The same procedure is then 
chosen, the integration between and e is performed 
using the series expansion and gives a contribution 



[h 2 t e 



a, a are not in agreement with our scaling analysis 
of PVT data on He 4 . Further investigations of the 
properties of the function (AIII.2), its range of validity 
and the sensitivity of the values of the parameters 
to variation in the form of the function are needed. 



+ & 4 



a + 2 



+ . 



(AII.6) 



whereas the integral between e and Xo is evaluated 
by standard numerical integration techniques. 

12. Appendix III 

Moldover's analysis of his He 4 specific heat data 
on 7 isochores is based on the assumptions that 
T(d 2 fJLldT 2 ) Pc is a constant (independent of T) and 
that the density independent part of the specific heat 
possesses a weak temperature dependence: 

-T[d 2 /jL(p c , t)ldt 2 ]-T(d 2 AoldT 2 ) = a + bt. (MILL) 

The singular part of the specific heat in the one- 
phase region, the term — T | Ap | _a//3 a"(x) in eq (2.21), 
is chosen as 

-T\Ap\- a !Va'\x)=A\&p\-«ie[l+x\x\ x - i ]-«l x . (AIII.2) 

The complete expressions for the specific heat that 
are used to fit the experimental data are as follows: 

pC v = Dkp + a + bt + A\Ap\-<*lV[\ + x\x\ N - l Y<*l N , 
one-phase region 



During the course of this research we have bene- 
fited from discussions with and preprints from Dr. R. B. 
Griffiths, Dr. M. E. Fisher and Dr. L. R. Wilcox. 
We made use of prepublication results kindly for- 
warded to us by Dr. M. R. Moldover, Dr. P. R. Roach 
and Dr. M. H. Edwards. 

We acknowledge the enthusiastic support of Dr. 
M. Klein and Dr. R. D. Mountain. 

Mr. J. Gallagher capably handled the computer 
plotting of several of the figures. Dr. M. Cooper 
and Dr. M. Barmatz provided constructive criticism. 

One of us (M.V.M.) has been partially supported 
by a NATO grant. 



pC v =Dkp + a' + bt+A'(-t)- a \ 

two-phase region 

Z>, a, a', a, a', b, N, A, and A' are constants to be 
determined by the fit. 

The assumption AIII.2, allows to analyze simul- 
taneously all C v isochores, which is a very attractive 
procedure. However as Moldover points out, the as- 
sumed functional dependence of x is not in complete 
agreement with the scaling ideas; the most serious 
problem (in our opinion) being that it does not have 
a power series expansion around x = 0. 

With eq (AIII.3) the specific heat data on all iso- 
chores scale nicely. The numerical value of the con- 
stant D[=—T(d 2 /jLldT 2 ) Pc ] determined by the fit is in 
agreement with our value (see III). However, the 
constants a and a' are found to be different (a' is 
positive and a is negative) which implies that although 
the singularity in C v is assumed to be described by 
the term |Ap|~ a// V(x) in the one phase region and 
by a term A' {— t)~ a ' in the two phase region, an 
additional jump in C v is present on the phase boundary. 

The values of a' derived by this analysis, a =0.16, 
is in disagreement with the value a' = reported in 
his thesis and also with the value a —0.05 that we 
find fitting the data in the two-phase region using 
(AIII.3) with 6 = 0. 

The additional jump in C v and the high value of 
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